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The Road to Chaos is Filled with 
Polynomial Curves 


Richard D. Neidinger and R. John Annen, III 


1. INTRODUCTION. The bifurcation diagram (Figure 1) is much like an aerial 
photograph — there are many interesting things revealed but it’s hard to see the 
detail or to interpret the features. By superimposing a family of polynomial curves 
(Figure 2), this picture turns into a road map — marking trails, boundaries, mile- 
stones, and labeling points of interest. Both the curves and the diagram are 
generated by iterating the function x 2 + c as the parameter ranges across a 
horizontal c-axis from 0 to - 2. The initial iterates form the curves with the nth 
curve being the nth iterate of the critical value zero as a function of c. The 
bifurcation diagram (also known as the orbit diagram, the Feigenbaum diagram, or 
the road to chaos) is a plot of later iterates, showing the eventual state of the 
system. 



-2 — 1.5 -1 - 0.5 0 

Parameter c 


Figure 1. The bifurcation diagram of f c (x ) = x 2 + c. 

In the first half of the paper, curve intersections provide insight into the 
dynamics of the system. The roots of the curves are shown to mark periodic 
windows in the bifurcation diagram, usually of invisible width, indicating a remark- 
able preponderance of attracting periodic phenomena that contrasts with the 
seeming chaos in the diagram. Other intersections mark parameter values where 
the iteration is chaotic. A brief detour into the complex plane yields a simple 
algorithm that color codes the Mandelbrot set to indicate the periods of the 
attracting phenomena. The second half of the paper studies how curve inequalities 
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-2 - 1.5 -1 - 0.5 0 


Figure 2. The first six Q-c\irves reveal dynamics and shapes within the diagram. 


describe the shapes seen in the bifurcation diagram. We show why the high-num- 
bered iterates condense along the polynomial curves. The curves also form en- 
velopes that bound embedded, self-similar copies of the bifurcation diagram. 
Curve intersections specify the exact beginning and ending parameters for every 
envelope found in the diagram and, thus, mark where each window ends. Finally, 
Sarkovskn’s ordering can be seen in the order in which the curves first cross the 
axis. The ordering is argued without using Sarkovskil’s theorem. 

Let’s begin by defining terms and notation. For the standard dynamical system 
f c (x) = x 2 + c, the iterates are defined by x j+1 = f c (x), or equivalently x t = f‘(x 0 ), 
for some fixed x 0 . A sequence of such iterates is usually called an orbit and the 
dynamics of the system refers to the behavior of orbits. The bifurcation diagram, in 
Figure 1, displays the dynamics for different values along the horizontal c-axis by 
vertically plotting iterates, numbered 300 to 900, in the orbit of zero. 

Define the nth polynomial Q n (c) = /"(! 0) to be the nth iterate of f c starting 
from zero. Specifically, Qfc) = c, Q 2 (c) = c 2 + c, and Q 3 (c) = (c 2 + c) 2 + c. In 
general, Q n+l (c) = (Q„(c)) 2 + c. Thus, each Q n (c) is a 2' ,_1 degree polynomial. 
The graph of any Q n as a function of c will be called a Q-curve. The first six 
(3-curves are displayed in Figure 2 and portions of these curves are visible as a 
relative density of points in Figure 1. 

Q-curves are iterates of the critical value (in [PJS], p. 633, they are called 
critical-value lines). The important role of the critical value is found in Fatou’s 
Theorem: Every attracting cycle for a polynomial attracts at least one critical point 
[Br], Even if the iteration for a bifurcation diagram starts from a different x 0 , the 
eventual state reveals the same 2-curves. Indeed, the bifurcation diagram in the 
background of Figure 2 uses x 0 = 0.2, while the curves are iterates of the critical 
value zero. Typically, even bifurcation diagrams start with zero in order to avoid 
isolated values where x 0 may be a repelling periodic point. 

2. CURVE INTERSECTIONS REVEAL SYSTEM DYNAMICS 

2.1 Roots and windows. In the bifurcation diagram, the gray haze of iterates 
contains fascinating gaps, or vertical strips of white space, usually called windows. 
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The windows are formed whenever the iteration abruptly changes from a wide 
spread of values to periodic oscillation. Each window can be associated with the 
lowest period of oscillation that occurs within that parameter range. A few 
windows are visible in Figure 1 but the (2-curves mark many windows that are too 
narrow to be visible. 

A quick observation from Figure 2 is that the curve Q 3 crosses the axis in the 
region usually called the window of period 3. Likewise, Q s and Q b cross in 
windows of period 5 and 6, respectively. The natural conjecture is true: whenever a 
(2-curve meets the axis, it marks a periodic window in the diagram (with period of 
the lowest numbered (2-curve that crosses there) and all windows are marked in 
this way. Actually, each root is a parameter value with a superattracting periodic 
cycle. Recall that p is a periodic point of period n if /"( p) = p but flip) # p for 
0 < j < n. This p is repelling if K/")'(p)l > 1, attracting if |(/ c '')'(p)l < 1, and 
superattracting if K/")'(p)l = 0. 

Superattracting Root Theorem. Let n e N. The parameter r satisfies Q n (r ) = 0 and 
Qjir) # 0 for 0 <j < n if and only if iteration of f fix ) has a superattracting periodic 
point of period n. 

Proof: Fix r. By definition, Q„(r) = /"( 0). Thus, Q n (r) = 0 and Qfr) + 0 for 
0 < j < n if and only if zero is a periodic point of period n. In general, a periodic 
point is superattracting if and only if a critical point belongs to the cycle. Indeed, 
by the chain rule, (f n )'(x 0 ) = n-’Jo/'C*,) and the product is 0 if and only if 
/'(*,) = 0 for some i [Br]. Since zero is the only critical point of f r , the theorem 
follows. ■ 

Such parameter values do mark windows. Indeed, continuity implies that every 
parameter with superattracting period n is contained in an open interval of 
parameter values with attracting period n. Conversely, any maximal interval of 
attracting period n will contain a superattracting value. In fact, more is true: the 
derivative of the nth composition at a cycle point ranges from 1 to - 1 across an 
interval of attracting period n. This statement follows from the more general 
complex result of Douady, Hubbard, and Sullivan ([Br], p. 85) and can be observed 
in the diagrams of Devaney [D2]. A window of period n in the bifurcation diagram 
is an interval of attracting period n that is augmented on the left by a region of 
bounded iterates that we identify as an envelope in Section 3.2. 

The roots of the (2-curves indicate many windows that are not visible in the 
bifurcation diagram (Figure 1) or even in zooms on the diagram [D2]. In Figure 3, 
the graph of Q n alone marks a tremendous number of windows of period 11. 
Curve (2 n has no roots (except 0) to the right of -1.5, but the zoom shows how 
roots, and therefore windows, are concentrated near - 2. 

2.2 Roots and Mandelbrot buds. The Superattracting Root Theorem can be used 
in the complex c-plane to mark components of period n. This is analogous to 
marking windows on the c-axis. The polynomials <2„(c) = / c "( 0) make sense for 
complex c-values, even though the term “curve” is not appropriate in this context. 
At complex roots where these polynomials are zero, the system has a superattract- 
ing periodic point of period n. Around each root, the region of attracting period n 
is a hyperbolic component of the Mandelbrot set [Br], and appears as a bud 
sprouting off a larger component or as a cardioid of a (mini-)Mandelbrot set. We 
can color neighborhoods of the roots, where |/"(0)| < e, in order to show the 
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-2 -1.99 -1.98 -1.97 -1.96 -1.95 

Figure 3. Every root of Q n marks a period 11 window. 


period of the corresponding components. 

A program for plotting the Mandelbrot set is easily adapted to perform the 
color-coding of components. For each pixel c, the typical Mandelbrot set algorithm 
starts with z 0 = 0 and iterates f c until |z„| > 2, so that c is outside the Mandel- 
brot set, or n exceeds some maximum iteration, so that c is assumed to be in the 
Mandelbrot set. We simply add a third bail-out option: whenever |z„| < e, and 
n <. the number of colors to be used in period-labeling, stop and color the pixel 
with color n. Choose solid colors, for the original two stopping tests (as opposed to 
the usual practice of coloring according to n when |z„| > 2). A fairly large e, such 
as 0.2, is recommended to color “disks” that very roughly cover the buds that are 
being labelled. More importantly, it will label (the period of points in the main 
cardioid of) mini-Mandelbrot sets near the boundary that would not be visible and 
could even be smaller than a pixel. 

This might be called a quick and dirty algorithm. It’s quick because it’s simple to 
program and actually provides an early exit for many pixels that would otherwise 
reach the maximum iteration. It’s dirty because it’s difficult to say any thin g 
rigorous about colored pixels for arbitrarily or aesthetically chosen e’s. It is usually 
easy to identify the component that goes with a colored disk, or vice versa, but 
colors can extend into adjacent components. Colored disks can be smaller or larger 
than the corresponding components, even for the same n and e. Other algorithms, 
such as direct period detection from iteration or the spider algorithm, are more 
rigorous. 

Two examples of the beautiful and intriguing color figures that are pro- 
duced by this algorithm can be viewed through the World Wide Web at 
http://www.davidson.edu/academic/math/neidinger/road.html. Fascinating 
patterns appear in the sequence of periods found in the geometric patterns. Of 
course, zooms near the boundary are as entertaining as ever. A period-n bud, off 
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the main cardioid, is crowned by a decoration consisting of n spokes and each 
spoke contains a sequence of periods (covering cardioids of mini-Mandelbrot sets) 
from an equivalence class mod n . 

2.3 The intersection dichotomy. We now leave the complex detour and return to 
the road map where every intersection is interesting. There are only two types of 
intersections between Q-curves and they imply different system dynamics for the 
corresponding parameter values. Figure 4 shows many intersections of the first 
eight Q-curves on [-2, - 1.8] and labels one of each type. A tangency occurs at a 
superattracting root r, and a non-tangent intersection occurs at a Misiurewicz 
point s, which has chaotic dynamics. While the implications of a root are relatively 
straight-forward, the implications of a Misiurewicz point rely on references in the 
literature. 



2 Li . . « i .. ... l i * i i * * * i , l i i , d 

-2 r - li9 s -1.8 


Figure 4. The two types of intersections: all curves tangent or none tangent, with periodic or chaotic 
dynamics, respectively. 


A Misiurewicz point is a c-value where the orbit of the critical value zero is 
strictly pre-periodic, i.e., it eventually falls onto a periodic orbit not including zero. 
The proof that the corresponding dynamics is chaotic is very involved and uses the 
definition of chaos from ergodic theory (see references in the proof below). 
However, the follpwing intuitive argument, from [CE], p. 34, yields some insight 
into the chaotic behavior: At a Misiurewicz point, the eventual periodic orbit of 
zero is repelling and, By Fatou’s Theorem, no periodic attractor can exist. We now 
argue sensitive dependence on initial conditions, meaning that nearby points are 
eventually “pulled apart” through the iteration. This condition is hardest to fulfill 
for points near x = 0, since the derivative of f c (x) at zero is zero, and, hence, f c 
contracts points near zero. But eventually zero itself is mapped onto a repelling 
periodic orbit. Thus, the orbits of points near zero may initially move closer (to the 
orbit of zero) but will eventually come near the repelling periodic orbit and be 
pushed away. 

Intersection Dichotomy. At any intersection of two Q-curves, the c-value is either : 

(1) a root of some Q n ( let n be the smallest such positive number ), where every 
intersection of Q-curves is a tangency at one of n distinct points ( in equivalence 
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classes modn), and where the system f c has a superattracting periodic point of 
period n, or 

(2) a Misiurewicz point , where no intersection of Q-curves is a tangency, and where 
the system f c is chaotic. 

Proof: Let c be the value at the intersection of two Q-curves. By definition of Q i9 
the orbit of zero is Qfc). If c is not the root of any Q-curve, then the orbit never 
returns to zero but the intersection implies that the orbit is eventually periodic, so 
c is a Misiurewicz point. Thus, c is either a root or a Misiurewicz point. 

The claims about system dynamics are found, in case (1), in the previous 
Superattracting Root Theorem and, in case (2), in results of M. Misiurewicz found 
in [CE] (p. 155 and supporting text). There, it is shown that if c is a Misiurewicz 
point, then f c has an ergodic, invariant, probability measure that is absolutely 
continuous with respect to Lebesgue measure. This measure-theoretic condition is 
one way to define a chaotic dynamical system. The proof in [CE] applies to a class 
of functions, including g = f c /c for a Misiurewicz point c < —0.5. Ergodic theory 
techniques are used to construct an invariant measure v (meaning v(E) = 
v(g~ l (E)) for every measurable set E) which is absolutely continuous. Such a 
measure is shown to be unique and, hence, ergodic. Readers unfamiliar with this 
characterization of chaos might want to consider how an attracting periodic orbit 
would prohibit the existence of an absolutely continuous invariant measure. 

To prove the claims about tangency, observe that: an intersection and its 
tangency or non-tangency can be propagated to the next iterate (and hence all 
subsequent iterates) by the recurrence relations Q /+1 (c) = (Qfc)) 2 + c and 
Q!+i(c) = 2Qfc)Q\(c) + 1. 

First, consider a root c = r where Q„(r) = 0 and no smaller numbered curve is 
zero at r. Then, Q, l+1 (r) = 0 2 + r = Q x (r) and Q,' + 1 (r) = 2 Q n (r)Q' n (r) + 1 = 

1 = QiW- By the recurrence relations, Q n+i (r) = Qfr) and Q' n+i (r) = Q/(r) for 
all i > 1. So all these intersections are tangencies at one of the n points 
Qi(r), . • • , Q, ,(/*). There are no other intersections of Q-curves at r since these n 
points are distinct (otherwise, forward propagation would yield Qj(r) = Q n (r) = 0 
for some j < n which contradicts the assumption on n). 

Now, assume that c = s is a Misiurewicz point. Then, since zero is not itself 
periodic, s is not the root of any Q-curve. Let k be the smallest number such that 
Q k intersects another Q-curve at s, and let n be the smallest number such that 
Q k (s) = Q k+n (s). This is not a tangency, by the following argument. By the 
recurrence relation, {Q k .f<s)) 2 = (Q k+n - iW) 2 and, hence, Q k - X {s) - 
-Q k + n - x {s). In particular, k > 2 (otherwise, the argument of the previous sen- 
tence would yield 0‘ = Q k+n -i(s), which contradicts no root at s). Now, define 
G(c) = Q k+n _ 1 (c) + Q k - X (c). In the next paragraph, we will show that G'(s) & 0. 
Then 

QUnO) ~ Qk(s) = 2 Qk+n-l( s )Qk+n-l( s ) ~ 2G*-l(O0UO) 

= -2G*_ 1 (5)G'(j)^0. 

We conclude that Q’ k+n (s ) ¥= Q’ k (s), so the intersection is not a tangency. By 
the recurrence relations, G„ +i (^) = <2,0) and Q' n+i (s ) ¥= Q!(s) for all i > k. 
(We can propagate the non-tangency since 5 is not the root of any Q-curve.) Thus, 
all these intersections are non-tangencies at one of the n distinct points 
Q k (s), ..., <2*+n_i( s ). There are no other intersections of Q-curves at s by the 
assumptions of “smallest” for k and n. 


1996] 


THE ROAD TO CHAOS IS FILLED WITH POLYNOMIAL CURVES 


645 



The following algebraic argument, from p. 333 of [DH], shows that G'(s ) # 0: 
By the recurrence relation, G'(s)/2 = Qk+n-i^Qk+n-i^ + Qk-i^Qk-i^ + !• 
To show that this expression cannot be zero, the rough idea is to show that the 
combination of Q and Q ' values has a positive power of 2 as a type of factor and, 
thus, cannot be - 1. Specifically, each rational number has a 2-adic valuation given 
by the power of 2 in the natural factorization of the fraction. The field of rationals 
and the valuation can be extended to include the number s. Let A be the set of 
elements of the extension field with non-negative valuation and let m be the set of 
elements with positive valuation. Then, m is a maximal ideal in the ring A. Using 
obvious valuations, such as 1 £ m and 2^4 c m, we work with arithmetic mod m. 
Since s is the root of a monic polynomial with integer coefficients (i.e., an 
algebraic integer), s and all integer coefficient polynomials of s (including all Q 
and Q' values) are in A. In fact, every Q[(s) = 1 mod m, since Q[(s) = 
2Q g - 1 (s)Q’ i _ 1 (s) + 1. Thus, 

Qk+n-2(s)Qic+n - 2 ( s ) + Qk-2 ( 0 Qk-2 ( 0 = Qk+n- lO) + Q k _ 2 (s)modm. 

A property of s, from the previous paragraph, is that Q k + n -i(s) + Q k ~\(s) = 0. 
By the recurrence, Q k+n - 2 (s) + Q 2 - 2 (s) = -2s e m. It follows that ( Q k+ „- 2 (s ) + 
Q k - 2 (s)) 2 and, hence, Q k +„- 2 (s) + Q k - 2 (s) are in m. By the preceding modular 
equivalence, Q k+n - 2 (s)Q' k+n _ 2 (s) + Q k - 2 (s)Q k _ 2 (s) is in m and cannot be -1 
(which has valuation 0). Thus, G'(s ) ¥= 0. ■ 

Now, one can literally see periodic and chaotic phenomena mixed throughout 
parameter ranges where the bifurcation diagram just shows a gray haze. Look at 
the curves to the left of the point marked r in Figure 4; there is superattracting 
periodic dynamics at every intersection with the axis, while there is chaotic 
dynamics at every non-tangent intersection of curves. Natural questions arise about 
the density and/or measure of the attracting periodic and chaotic parameters. The 
known answers are deep results in the field. Recently, Swiqtek showed that the 
union of the attracting periodic parameters is an open, dense set in the interval [S] 
and, yet, Jakobson’s theorem shows that the set of chaotic parameters has positive 
Lebesgue measure [R]! 

3. CURVE INEQUALITIES SHAPE THE DIAGRAM 

3.1 Why curves appear in the bifurcation diagram. In the bifurcation diagram, 
Figure 1, the first few 0-curves appear as pathways that are well-traveled by later 
iterates. Once the iterates spread across a vertical range (for c < the Feigenbaum 
point - 1.401155 . . . , [C], p. 141), they seem to cluster around curves, particularly 
on one side of each curve. This density of iterates is especially clear in the 
histograms of [D4], p. 127, and [PJS], p. 632. Why do the 0-curves, formed by 
initial iterates of zero, emerge from these later iterates, even when x 0 ¥> 0? 

One way to explain this phenomenon is to show how iterates in a relatively wide 
band around the axis must map into narrow bands around the first few 0-curves. 
These bands are shown in Figure 5, between solid 0-curves and corresponding 
dashed curves, for c in [ — 2, —1]. The dashed curves are given by P„(c) = f "(e) 
(for e = 0.15 in Figure 5), corresponding to 0„(c) = f"( 0). If any iterate x ( falls in 
the band of width 2e around the c-axis, then the next iterate will be in the band of 
width e 2 between Q 1 and P lt which condenses iterates above Q v Specifically, if 
0 < |x,| < e, then, by squaring and adding c to each side, 0/c) < x i+1 < P x (c). 
Assuming P x (c ) < 0 and again applying f c , yields 0 2 (c) > x i+2 > P 2 (c). Since 
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P 2 (c) > 0 for all c < —1.4, we can “propagate the bounds” to the next iterate, 
23(c) > x i+3 > P 3 (c); so that the density of points will occur below Q 3 . In the 
range of c-values where Q 3 (c) ^ 0, the next iterate jc i+4 will fall above Q 4 but, in 
the c-range where P 3 (c ) > 0, it falls below Q 4 . The bands should appear con- 
densed when the width I2;(c) — /^(c) is smaller than the initial 2e, which happens 
for low i and small e, though the width can grow as c approaches —2. We can 
estimate \Q i+1 - P i+1 | = \Qf - P 2 \ < (I2,| + IP; I) 1 2, - P, \ <, 4| 2, - P,|. By induc- 
tion, | Q i+1 — P i+ 1 < 4‘e 2 . This phenomenon emphasizes one side of low-order 
2-curves (lower should be darker); it can be observed in the bifurcation diagram in 
Figure 1. 



Figure 5. Iterates near the axis are mapped into concentrated bands around curves. 


A geometrical view of the dynamical system can explain why the whole region 
below Q 3 is darker than the region above it, as c ranges between about — 1.8 to 
— 1.4. The bifurcation diagram on [—2, — 1] can be divided into two regions, in the 
(c, x)-plane, depending on whether or not Ul < 22(c)- If (c, x) satisfies 0 < |x| < 
22(c), applying f c yields 2i(c) < f c (x) < Q 3 (c ); i.e., it maps the symmetric region 
by folding it over along the axis, stretching the axis down to 2i> and stretching 22 
down to Q 3 . The remaining portion of the bifurcation diagram satisfies 2i(c) < x 
< -2 2 (c), so that .applying f c yields 2 2 ( c ) ^ / C U) ^ 2 3 (c); flipping it over and 
filling in above 23- In the range between about -1.8 to —1.4, this motion 
condenses the first region and expands the second. 

3.2 Envelopes bound copies of the diagram. On the road map given by curves on 
top of the bifurcation diagram (Figure 2), 2-curves form boundary lines separating 
regions that contain iterates from those that don’t. In fact, these curves bound the 
shape and identify the location of all embedded, self-similar copies of the dia g ram 
In general, whenever a curve Q n crosses the axis, the curves Q n and Q 2 „ bound a 
self-similar copy of the bifurcation diagram consisting of every nth iterate. The 
other iterates fall in similar envelopes, also bounded by 2-curves. 

Figure 6 shows three such sets of envelopes. First, the entire diagram is 
bounded by 2i below and Q 2 above. Second, all the even iterates are bounded 
between Q 2 and Q 4 , from the root of 2 2 to the point where Q 4 meets —Q 2 . In 


1996] 


THE ROAD TO CHAOS IS FILLED WITH POLYNOMIAL CURVES 


647 



this parameter range, the odd iterates fall between Q 1 and Q 3 . Finally, a quarter 
of the iterates are bounded between Q 4 and Q s , from the root of Q 4 to the point 
where Q s meets — Q 4 . Four envelopes can be seen in this parameter range. 



-2 -1.5 -1 -0.5 0 


Figure 6. g-curves bound nested copies of the diagram. 

Figure 7 zooms in on the period three window, around the axis, to show an 
envelope formed when Q 3 crosses the axis. The envelope abruptly ends at the 
Misiurewicz point where -Q 3 meets Q 6 . Within this range, every third iterate falls 
between Q 3 and Q 6 , while the other iterates are constrained to the region between 
<2j and Q 4 or the region between Q 2 and Q s . This containment is being observed 
when people speak of a “period 3 window.” This is different from the interval of 
attracting period 3. Actually, if the envelope is on the interval (s, r ) and the 
interval of attracting period 3 is (a, b), then the window is the union (s, b). (The 
right endpoint, b, of the attracting interval is a point not described by (7-curve 
intersections.) Each of the tiny, nested windows seen in Figure 7 contains even 
smaller envelopes of Q-curves. 



-1.79 -1.78 -1.77 -1.76 -1.75 

Figure 7. A zoom into the period 3 window shows an envelope of g-curves. 
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The following theorem describes all such envelopes. The phrase x lies between y 
and z means y<x<z ox z<x<y. 

Envelope Theorem. If Q n (r) = 0 and Qj(r) ¥= 0 for 0 < j < n, let s = 
max{c < r: Q ln (c) = -Q n (c)}. IfcE: [s,r] and \x Q \ < \Q n (c)\ , then x kn+m lies 
between Q m and Q m+n for all m e {1, 2, . . . , n] and k e {0, 1, 2, . . . }. In particular , 
Qkn+m H es between Q m and Q m+n on [$, r]. The regions between these pairs of curves, 
Q m and Q m+n for each ra e { 1,2, , n), are disjoint from each other on (s, r]. For 
later application, we also note that, ifm¥^n, the region between Q m and Q m+n never 
meets the axis on [s, r\. 

To prove the Envelope Theorem, we “propagate the bounds” from one iterate 
to the next. However, this requires either that the iterate must fall between - Q n 
and + Q n or that both the upper and lower bounds must have the same sign. This 
will be reduced to the requirement that \Q ln \ < \Q n \, a fact that is fairly obvious in 
the graphs. This and other essential properties of the “horn shape” formed by the 
graphs of Q ln and Q n (see Figures 6 and 7) are summarised in the following 
lemma. The claims following “moreover” are used in the next section. 


Horn Lemma. If Q n (r ) = 0 and Qj(r) ¥= 0 for 0 < j < n, then there exists s < r 
such that Q 2n (s) = - Q n (s ) and \Q 2n \ < \Q n \ on ( s,r ). Moreover, there exists 
r' e (s, r) such that 0 2 «( r ') = 0 , Q ln Qn < 0 on [s, /) and Q ln Q n > 0 on (r', r). 


This lemma is deceptively innocent-looking — it embodies deep properties of the 
quadratic map that we will not prove. The “horn shape” is a real-parameter 
consequence of Douady and Hubbard’s mini-Mandelbrot sets for the complex 
parameter [DH]; they prove that only complete (not partial) copies of the set occur. 
The horn shape is also equivalent to a pattern observed in graphical analysis (used 
to produce web-diagrams), as described in the following paragraph. 

Consider the graph of y = /"(x) near the origin in the Ay-plane. If Q n (r) = 0, 
zero is a superattracting periodic point of period n for f r , so that y = f"(x) has an 
extreme value at the origin. In fact, for a range of c-values near r, the graph 
resembles a parabola in a small region around the origin. We examine the 
concave-upward case shown in Figure 8. As c decreases from r to a value that we 
call s, the vertex moves monotonically downward until it becomes as large, in 
absolute value, 'as the coordinates of the intersection with the diagonal in the first 
quadrant. This diagonal intersection is a repelling point that forms the corner of a 
square of trapped iterates, as in [Dl], p. 132. The properties of the Horn Lemma 
can now be deduced from this family of graphs. On each graph of y = /"(x), the 
y-intercept is Q n (c) = / c "( 0). Now, Q 2n (c) = f' c \Q n (c )) is visualized by taking 
the y-intercept horizontally to the diagonal and then vertically to the curve. In 
Figure 8, the resulting values of Q 2n (c ) are connected, creating a curve parameter- 
ized by x = Q„(c) and y = Q 2n (c) for c in [s, r]. This parametric curve displays all 
of the properties stated in the Horn Lemma. This is not a proof, but merely shows 
that the lemma is equivalent to the claim that the parabola-shape always moves 
monotonically until reaching the edge of the trapping box. Graphical analysis also 
shows the dynamical significance of the parameter s — beyond s, the parabola 
pokes out of the box and the containment (or envelope) of iterates ends. We now 
return to the c-axis for the proof of the containment between Q-curves. 
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Figure 8. Five graphs of y = /"(*) lead to the Horn Lemma as shown in the thick curve parameterized 
by x = Q„(c) and y = Q 2 „(c\ 


Proof of the Envelope Theorem: (Assuming the Horn Lemma.) Fix r and n such 
that Q n (r) = 0 and Qfr) #0 for 0 < j < n. By the Horn Lemma, let s < r 
such that Q 2n (s) = ~Qn ($) an d < I2„(c)l for c e (s, r). Clearly, s = 

max{c < r: Q ln (c ) = - Q„(c)}. 

Fix any c e [s, r] and consider the orbit from any |x 0 | < \Q„(c)\. Apply f c to 
0 < |x 0 | < \Q n (c)\ to get Qfc) <x 1 < Q ll+1 (c). The containment asserted by the 
theorem continues because of two propagation principles. (1) If x lies between Q m 
and Q m+n for some m e {1,2, . . . , n - 1}, then f c (x) lies between Q m+l and 
Qm+i+n- (2) If x lies between Q n and Q ln , then f c (x) lies between Q x and Q 1+n . 
The second principle uses the fact that I22n( c )l < \Q,fc)\ for c e [j, r]. If x lies 
between Q n and Q 2n , then 0 < |x| < |Q„(c)| and, again by applying f c , Qfc) < 
f c (x) < Q ll+1 (c). To prove the first principle, suppose x lies between Q m and 
Q m+n . One possibility (out of two) is that Q m (c) <x < Q m+n (c). In the following 
paragraph, we show that Q m (c ) and Q m+n (c) agree in sign. If they are both 
positive, then Q m+X (c) < f c (x) < Q m+1+n (c). If they are both negative, then 
Q m+1 (c)>f c (x)^Q m+1+n (c). In either case, f c (x) lies between Q m+1 and 
Qm+i+n • Th e other possibility is that the original containment was Q m+n (c) < x < 
Q m (c ), and the argument is identical. To complete the proof of containment, it 
remains to show that Q m and Q m+n are nonzero and agree in sign on [s,r]. We 
will show a bit more as a bonus. 

We claim that the curves Q v Q 2 , ■ ■ ■ , Q 2n -i never meet the axis in (s, r ) and 
that none of the curves Q v Q 2 , Q 2n intersect in (s, r). At the root r, the curves 
pair up, mod n, at the n distinct points (by the Intersection Dichotomy). Thus, the 
claim implies that the regions between these pairs of curves, Q m and Q m+n for 
each m e {1, 2, . . . , «}, are disjoint from each other on (s, r]. The claim also 
implies that, for each m e {1,2 ,...,« - 1}, the pair Q m and Q m+n are nonzero 
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and agree in sign throughout (s, r\ Since s is a Misiurewicz point, it can’t be a root 
of Q m or Q m+n - So, in fact, Q m and Q m+n are nonzero and agree in sign 
throughout [s, rj. We prove the claim by contradiction. Let B be the set of 
counterexamples, i.e., B = {c e ( 5 , r): <2,(c) = Qj(c) for some i < j <, 2 rt, or 0 = 
Qj(c ) for some j < 2 «}, and suppose B is not empty. B is finite, since every <2, is a 
distinct polynomial, so we can let b = max(B). Since Qfb) = Qj(b) or 0 = Qfb), 
propagate (apply f b ) until Q k (b) = Q ln (b ) for some k < In. Now, consider the 
continuous curves \Q k (c)\ and \Q n (c)\ on (b, r ). At b, \Q k (b)\ = \Q 2n (b)\ < \Q n (b)\. 
By this inequality, k*n. At r, \Q k (r)\ > 0 = \Q n (r)\. By the Intermediate Value 
Theorem, there exists some d in (b, r ), such that \Q k (d)\ = \Q n (d)\. Propagate this 
equality once to get Q k+l (d ) = Q n+l (d). Thus d e B, which contradicts the 
assumed maximality of b. ■ 

The Envelope Theorem’s hypothesis that |x 0 | < \Q n (c)\ is necessary to have 
containment for all c in [ 5 , r] but, in practice, is not needed for most c-values. The 
bifurcation diagram, in the background of Figure 7, was created with x 0 = 0 to 
satisfy this hypothesis, which essentially requires that the iteration start inside an 
envelope. For other x 0 values, the diagram looks almost the same but there can be 
isolated exceptions where iterates leak out of the envelopes. Outside the en- 
velopes, an x 0 value can be, for instance, in a repelling periodic orbit for various 
values of c (with different periods) in the range of the envelope. At such c values, 
the iteration never enters the envelopes; nearby, it enters only after many itera- 
tions. Still, for most c-values, the iteration is rapidly attracted down into the 
envelopes. In fact, for a given c in the range of an envelope, the set of those values 
of x 0 with bounded orbits that are not attracted into the envelopes is a set of 
Lebesque measure zero. 

3.3 Sarkovskii’s ordering. The Q n curves also reveal the pattern known as 
Sarkovskn’s ordering. This ordering of N is as follows: 

3>5>7> ••• >2-3>2-5>2-7> ••• >2 2 -3>2 2 -5> ••• >2 3 -3>2 3 -5> ••• 

> ••• > 2 4 > 2 3 > 2 2 > 2 > 1 . 

Sarkovskn’s theorem [D2] applies to any continuous F: IR -» R and states that if F 
has a periodic point of period n and if n> m, then F also has a periodic point of 
period m. Although this theorem applies to f c for a fixed parameter c, it’s not 
immediately clear how this relates to Q n across a range of c-values. We will not 
use Sarkovskn’s theorem but will establish a relationship by directly studying the 
(2-curves. 

Root Ordering Theorem. Define r n = max{r: Q n (r ) = 0 and Qfj ) # 0 for 
0 <j < n). Then r n < r m if and only ifn>m. 

Proof: Every r„ exists since, for n > 2, Q n (- 2) = 2 and Q n < Q 2 < 0 on (r 2 , r^. 
By bifurcation theory (or the Horn Lemma and the Envelope Theorem), 
••• < r 8 < r 4 < r 2 < r x and if i is not a power of 2, then r, is to the left of this list. 
For each power of two, p = 2‘, there is a corresponding s p < r p described by the 
Horn Lemma. Also, s p < s lp < r 2p < r p . These inequalities follow from the Horn 
Lemma; s p can’t be in (s 2p , r 2p ) since on this interval \Q 4p \ < \Q 2p \ and yet, at s p , it 
can be argued that Q 4p (s p ) = Q 2p (s p ). We conclude that s 1 <s 2 <s 4 <s 8 < 
••• < r 8 < r 4 < r 2 < r v To complete the proof, we now prove that, for each power 
of two p = 2 1 , s p < r 2p < r 5p < r lp < < s 2p . The reader can use Figure 9 to 

follow along for p = 1. 
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uses: . , i 1.1 1 ,, l , , i , , , , 

■*i s 2 r i r 2 r l 


Figure 9. Sarkovskil’s ordering appears as odd curves (shown in black) cross the axis between jj and s 2 . 


Let k e {3, 5, 7, . . . }. First, we have already argued that bifurcation theory (or 
the Horn Lemma and the Envelope Theorem) ensures that r kp < r 2p . On [s 2p , r 2p \, 
the Envelope Theorem for n = 2p says that Q kp lies between Q p and Q 3p and the 
region never meets the axis. Thus, r kp < s 2p . We now establish that the roots do 
occur in the interval ( s p , s 2p ). We will repeatedly use the fact that, at s 2p , the sign 
of each Q kp agrees with the sign of Q p , a fact implied by the observed containment 
of Q kp curves. 

By the Horn Lemma for n = p, Q p and Q 2p have opposite signs throughout 
[s p , r 2p ), in particular throughout [s p , s 2p \. At s p , -Q p (s p ) = Q 2p (s p ) and iteration 
yields Q 2p (s p ) = 0 3p (s p ). Thus, the sign of 0 3p changes from the sign of Q 2p at s p 
to the sign of Q p at s 2p . By the Intermediate Value Theorem, there exists a largest 
root r 3p of 0 3p in (s p , s 2p ). Now, at r 3p , 0 5p (r 3p ) = Q 2p (r 3p ). The sign of Q Sp 
changes from the sign of Q 2p at r 3p to the sign of Q p at s 2p . By the Intermediate 
Value Theorem, there exists a largest root r 5p of Q Sp in (r 3p , s 2p ). Again, at r 5p , 
Q lp (r 5p ) = Q 2p (r 5p ), so that the sign of Q lp changes on (r 5p , s 2p ). By induction, 
s„ < r 3 P < r s P < r lp < - < s 2p . ■ 

4. CONCLUSION. Many of the most intriguing features along the road to chaos 
are clarified by looking at the road map showing g-curves. For years, people have 
studied the wonders of the bifurcation diagram: there appears to be a region of 
chaos that is interrupted by windows of periodic attraction and that has shadowy 
curves from varying density of iterates; the entire diagram appears in smaller, 
similar copies nested within the diagram. Each of these phenomena is better 
understood if one also studies the Q-curves. The roots of Q-curves show that 
periodic windows are much more prevalent than the diagram suggests. On the 
other hand, the non-tangent intersections of Q-curves point out at least some 
points where one can confidently say that the system is chaotic. The Q-curves are 
the curves that appear in the diagram, both as the shadowy curves of higher 
density and as the boundaries of the diagram and the nested copies. Indeed, any 
root of a 0-curve begins, and a corresponding intersection ends, a nested copy of 
the diagram. Even Sarkovskifs ordering and a color-coded (according to period) 
Mandelbrot set arise from this study of 0-curves. 
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Most of our proofs use only elementary mathematics and most of our arguments 
remain on the parameter axis, rather than switch to web diagrams on the graph of 
the quadratic. The key ideas are very accessible — periodic attraction at roots uses 
the chain rule, tangencies and envelopes use elementary algebra on equalities and 
inequalities, and Sarkovskn’s ordering uses the Intermediate Value Theorem. 
However, the properties of Misiurewicz points and the Horn Lemma do rely on 
more advanced work in ergodic theory, algebra, and complex analysis. There is 
ample graphical evidence for the Horn Lemma, which plays a crucial role in the 
Envelope Theorem and Sarkovskifs ordering. 
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